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ABSTRACT 


This review article critically outlines and discusses the main issues of 34 methods which have been 
developed and validated over the past 35 years in order to determine with an acceptable accuracy and 
reliability fundamental parameters of solar cells. This review covers methodologies which deal with 
current-voltage characteristic (I-V) analysis either theoretically through elaborated models and/or 
treated graphically. Methodologies based on the theoretical analysis of the I-V characteristics using 
the one or two diode model are discussed. The investigation on the I-V characteristics is processed via 
statistical functions, non-linear regression and stochastic models. A second family of methods to 
determine the solar cell electric parameters comprises the ones which deal with the graphical treatment 
and analysis of the I-V characteristics which are measured at different environmental conditions. To the 
third family belong the methods which use a mix approach of theoretical analysis of the I-V 
characteristics through modeling on one hand and the graphical analysis of their experimental 
configuration, on the other. The paper comments on each of the 34 methods and provides pros and 
cons for the determination of the fundamental electric parameters of solar cells. 

© 2013 Elsevier Ltd. All rights reserved. 


1. Introduction 


The determination and knowledge of the influence of solar cells 
parameters upon their efficiency become crucial for the optimiza- 


The main targets in the energy policy and the research 
priorities and efforts drive towards increasing the efficiency of 
the renewable energy systems available, and to the design of cost 
effective technologies. 

The production of photovoltaic energy has undergone a sustainable 
growth rate, with maximum points in the three oil crisis, when the 
need of alternative types of energy sounds prominent. The objective 
established by the European Union is that until 2020 the photovoltaic 
energy represents 12% of the total electric energy production [1]. 
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tion of fabrication processes and for the scientific research [2]. The 
possibility to determine the important parameters of the solar 
cells plays a major role in evaluating the solar cells, in controlling 
their quality, and also in the fabrication and power performance of 
reliable solar panels. 

The most important parameters of solar cells can be deter- 
mined by using the current-voltage (I-V) characteristic which is 
shown in Fig. 1 and by analyzing their equivalent circuit [2]. These 
parameters are: lpn - the photogenerated current, Isc — the short 
circuit current, Voc — the open circuit voltage, n - the ideality factor 
of diode, R, — the series resistance, R,, — the shunt resistance, 
I, - the reverse saturation current, P,, — the maximum power, Vm - 
the voltage at Pm, Im — the current at P,,, FF — the fill factor and 
n — the efficiency. 
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Fig. 1. The I-V characteristic for solar cell. 


Fig. 2. The equivalent circuit for a solar cell - the one diode model. 


The fill factor and the efficiency of a solar cell can be calculated 
by the expressions 


Pm _ ImVm 


FF= = 1 
IscVoc IscVoc ( ) 
Pm 
n= Ary (2) 


where A represents the area of solar cell and I; the solar irradiance 
with spectrum AM 1.5 and T=25°C, that is Standard Test 
Conditions (STC). 

The solar cell electric equivalent circuits and the corresponding 
mathematical models take into consideration the number of 
diodes in the circuit. Fig. 2 presents the equivalent circuit of the 
one diode model which is mostly utilized to characterize the 
solar cells. 

The mathematical model which describes the one diode model 
is characterized by the following equation [3]: 


V+IR; 
Rsh 


ael Ga 1) (3) 
where Vy represents the thermal voltage, V;=kT/q, k denotes the 
Boltzmann constant, T is the temperature of solar cell and q is the 
electronic charge. 

In case the electronic conduction mechanisms within a solar 
cell are considered separately, this can be modeled by a power 
source and a diode in parallel for each mechanism in turn: 
diffusion, generation-recombination and thermionic. In this case 
the above mathematical model is modified by adding an expo- 
nential term for each mechanism. Thus, the two diode model, 
which considers the diffuse and recombination mechanisms, is 
described by Eq. (4) [4], while Eq. (5) [4] describes the three diode 
model: 


+IRs V+ IRs. 
T= (e1) io (e mV; -1) aie (4) 
S| 


where d and r indices represent the diffusion and the generated- 
recombination mechanisms, 


V+IRs V+IRs IRs 
T= Iph —loa (ert) —lor (e "Vr 1) Tot (err 1) Vat (5) 
Rsh 


where t index denotes the thermionic mechanism. 


In general for the silicon solar cells it is useful to take into 
consideration the two diode model, but in case of thin film cells 
(heterojunctions) the second exponential term has a small influ- 
ence, which can thus be neglected [5]. 


2. Methods used 


The measurement of the I-V characteristic of the solar cells can 
be done in the laboratory or in natural sunlight conditions. In the 
former case the measurements can be performed under dark or 
illuminated conditions, in the latter case solar simulators being 
used [6-25]. 

Most of the methods used to determine the parameters of the 
solar cells use the one diode model, as the interpretation of the 
equation describing the mathematical model is simpler. 
The methods to determine the parameters aforementioned are 
briefly presented below. 

Method 1. Chan et al. developed the analytical five point 
method [9]. This method uses the one diode model for solar cells 
and allows the determination of the following parameters: Ipn, Io, 
n, R; and Rsn (the five parameters - P) under illuminated condi- 
tions using the values of Iseo Voc, Im, Vm, the slope at open-circuit 
point R,, and the slope at short-circuit point R,,. which are 
obtained from the I-V characteristic. The values of Rsno given by 
Eq. (6) and R,o by Eq. (7), can be calculated by linear fit of I-V 
characteristic around the short circuit current point and around 
the open circuit voltage point respectively. The shunt resistance is 
calculated using the following expression: 


dV 
Reh = Rho al (T), oe (6) 
dV 
ss (ar) E D 
The ideality factor of the diode can be calculated by 
A 

Aee i 
where A, B and C are calculated as follows: 

A= Vm + Rsolm—Voc (9) 

Vin Voc 
B= In (Is Re, In) In (1s re) (10) 
C= a (11) 
Isc— Ree 


The other parameters lo, Rs and Ip, are obtained using the follow- 
ing equations: 


Voc Voc 
= (i = exp( a) tal 
nV; Voc 
Rs = Rso T, exp( wie) (13) 
_ Rs IscRs 
Inn =lee( 1+) +1o(exp (2) -1) (14) 


Method 2. Ishibashi et al. proposed a new method for the 
derivation of all solar cell parameters. The parameters are obtained 
by a single I-V characteristic measured under illumination condi- 
tions, using the one diode model with only one assumption, 
Eq. (15), which is true for the most types of solar cells [10]. 


a= exp (“os ) a1 (15) 
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Using the equations obtained from Eq. (3) for the coordinates (0, 
Ic) and (Və, 0), and the assumption Eq. (15) the following 
expression is obtained: 


dv nvr 4 
dl Isc—I—[V—Rs(Isc—D—nV7]/Ren 


Rs (16) 


Substituting I=/,. and V=O into Eq. (16) the relation becomes 


dV , 
-(T), TP+ Ra a7 


where Rs<Rsn. Ishibashi et al. introduced two temporary para- 
meters in Eq. (16), Rso and n, to determine the series resistance 
and the ideality factor of diode. These two terms are used in an 
iterative process [7]. By plotting — dV/dI vs. {Isce—I-[V —Ros(Ise—D— 
NoVr]/Rsn}—!, n can be obtained as a slope and R, as a y-intercept. 
The iterative process begins with R,,.=0 and ny=0 and ends when 
|5|=0, where 6 is given by the following formula: 


(Rs—Rso)Usc—D +(n—No)Vr 
V—Rso(Isc—D—NoVr 


b= (18) 


Therefore, the other parameters, I, and Ip, can be obtained 
using the following expressions: 


Voc—IscRs Voc 
m= (1« Rsh ) exp( ws) (19) 
= Voc Voc 
Iph =R, th (exp (že) -1) (20) 


Method 3. Chegaar et al. developed a simple method to extract 
the following parameters of the solar cell: Rsm, Rs, Io and n. The 
parameters are obtained under illuminated conditions using a 
single I-V characteristic, the one diode model, an auxiliary func- 
tion F(V) and a fitting routine [2]. The shunt resistance is obtained 
using Eq. (6). An auxiliary function is introduced by Chegaar et al. 
to determine the other three parameters: 


F(V) = V—Va In(Ipn—1) (21) 


where V, represents an arbitrary value of the voltage. Using Eq. (3) 


without the last term and I instead of V, Eq. (21) becomes 
where a=—R;s, b=nVr and c=—nVy7I, 


(22) 


FM = al+b In(Ipn—I) + Co 


The series resistance, the ideality factor of the diode and the 
reverse saturation current can be determined by fitting the curve 
F(T) vs. I. The curve is obtained by calculating F(I) for each point (V, 
I) at a constant temperature and for a value of Va. Chegaar et al. 
take into consideration integer values for V, between 1 and 5, [2]. 
The photogenerated current is considered approximately equal 
with the short circuit current. 

Method 4. The parameters of the solar cell such as Rsn, Rs, lo and 
n can be calculated with the method developed by Tivanov et al. 
[11]. Considering the mathematical equation for the one diode 
model, which is represented by Eq. (3), approximating Ipn ~ Isc and 
analyzing the I-V characteristic of solar cell at a fixed level of 
illumination, its parameters can be determined by the following 
equations: 


97 1/2 
R= : (a-b) 4 Pa b)+ ee) (a+b) +2 (23) 
2 Tee Tse Isc 


Voc 
Ra =- (24) 
O ptk arktis 


Isc Voc Voc 
= = ex 
Rye Pny : (25) 
di\~! dI\ ~! 
= yt+2_ = — 
a ai a= (5), o b= (T) a ai 


The assigned particular value of n is verified by comparison with 
the value obtained with the following equation: 


(Vm—ImRs)(Isc+lo—Im—Vm/Rsn) 


= 27 
Vr(Im—Vm/Rsn) € ) 
In the particular case when a> —oo Eqs. (23)-(25) become 
nV I 
R; = |b] T L R>, lo= x (28) 
sc Y 


Method 5. Easwarakhanthan et al. developed the vertical 
optimization method. The one diode model was used for solar 
cell. This method allows the determination of the five parameters 
of solar cells using the non-linear least squares optimization 
algorithm [12]. The modified Newton method with the Levenberg 
parameter is used for the algorithm. The objective function S in 
case of the study of a solar cell is 


N 
S(P)= > (li—I(V;, P))? (29) 
i=1 


where l; V; are values on the I-V characteristic whose mathema- 
tical model is given by Eq. (3) and P represents the set of the five 
parameters. The goal is to minimize the S function. 

Method 6. The electric conductance optimization method was 
developed by Chegaar et al. [13]. This method allows the determi- 
nation of the five parameters of a solar cell using the non-linear 
least squares optimization algorithm considering the one diode 
model for the solar cell. The method is based on the modified 
Newton method with the Levenberg parameter used. This method 
considers the I-V characteristic and the slope of this curve. The 
objective function S; in case of the solar cell is 


N 

Si(P)= $ (G-Gi(V; Ii, P)? (30) 
i=1 

where G is the conductance calculated with Eq. (31), I; and V; are 

the values on the I-V characteristic, G; represents the computed 

conductance and P represents the set of the five parameters. The 

goal is to minimize the S4 function. 


Ipn—lo—(V + IRs Gch +1V 7Goh 
nV; +R; (Ipn—Lo—(V + IRs )Gsh + NVTGsn) 


where Gn = E 
sh 


(31) 


Method 7. Chegaar et al. proposed the simple conductance 
technique based on Werner’s method | 12,13]. The five parameters 
can be determined using the one diode model and Eq. (3) in two 
steps. In the first step the shunt conductance Gsp is determined by 
simple linear fitting of the reverse bias characteristic for a large 
negative bias voltage. The second step consists of rewriting Eq. (3) 
in forward conditions for V+R,J>kT and calculating the conduc- 
tance G=dl/dV: 


G_41R 
Ipn—I g “nV; nVr 


(32) 


Using the plot of Eq. (32), the ideality factor n of the diode can be 
determined by the y-intercept and the series resistance R, by the slope 
respectively. For the photogenerated current, it is highly acceptable to 
approximate lpn ~ lsc. Taking into consideration the I-V data corrected 
using the effect of the series resistance, I, can be determined by the 
standard method based on the forward I-V data [4]. 

Method 8. Haouari-Merbaha et al. developed a method to 
determine the parameters of solar cells based on the two diode 
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model [15]. The I-V characteristic of a solar cell measured under 
illumination conditions is split in two regions by the resistive line 
with slope V,,/I;-, the first region is near the short circuit current 
and the second region is near the open circuit voltage. These two 
regions of the I-V characteristic are fitted as follows: the first 
region as a function of the current and the second as a function of 
the voltage. For the first region Eq. (4) can be modeled as Eq. (33) 
and for the second as Eq. (34): 


I= Ise—G, V—a3 (I, V) (33) 


V = Voc—IR? —b3h° (I, V) (34) 


where the superscript “s” is for the first region, the superscript “o” 
is for the second region, h'(I,V) represents the curvature of the I-V 
characteristic, i can be either s or o [14]. The coefficients (Isc, Gsh?) 
can be obtained by fitting the first region near short circuit, with 
Eq. (33) and the coefficients (Voc, R°) can be obtained by fitting 
the second region near open circuit, with Eq. (34). The set of 
parameters (Ipn, Gsn= 1/Rsn, Tod, lor, Rs) can be determined by solving 
the system with five non-linear equations obtained for: Ise Gsn°, 
Voc, Rs° and the curve passes through the boundary point, which is 
around the maximum power point [15]. The ideality factors of the 
diode take the theoretical values ng=1 and n,=2. Solving the 
system starts with the set of parameters equal with zero and the 
iterative process is done when the parameters converge. 

Method 9. The method proposed by Ortiz-Conde et al. to 
determine the parameters of a solar cell is based on the Co- 
content function CC | 16]. The one diode model and Eq. (3) are used 
to characterize the solar cell. The I-V characteristic is obtained 
under illumination conditions. Using the Lambert W function and 
the Co-content function CC, Eq. (3) can be rewritten as 


CCA, V) = Cvi V + Cn Isc) + Crvi VU Isc) + Cy2V? + Cia (I-lsc)* 
(35) 
The equation coefficients Cy, Cn, Cvz2 and Cp can be determined 


through fitting Eq. (35) using the CC function. The CC function is 
numerically calculated from the measured data as follows: 


V 
Ch V) = { (I-Ige)dV (36) 
0 


Provided the equation coefficients are determined, the parameters 
of the solar cell can be calculated by 


1 
Ren =e, (37) 
14+ 16C Cp 1 
R= 4v2 (38) 
Cji ( T+ T6Cy2C2—1) +4CnCv2 
iz (39) 
AV1Cy2 
(Cv +Ise)( 14 16Cv2Ce +1) 
ph = 5) 2Cr Cv2 (40) 
I—(V-IRs)/Rsn + Iph an 


° exp((V—IRs)/nV7)—1 


Method 10. Zhang et al. developed a method to determine the 
parameters of the solar cell based on the one diode model, the 
single J-V characteristic under the constant illumination level and 
the Lambert W function |17]. The implicit Eq. (3) can be trans- 
formed into an explicit equation with five unknown parameters 
Ipm lo, n, Rs and Rsn by using the Lambert W function. Using Eq. (3) 
in the particular points (0, I) and (V,,, 0) and the assumption 


given by Eq. (15), a simplified explicit equation can be obtained: 


f= "g lambert W Ea (Isc K sa) exp CAR (Rslsc + nik) | 


Vy RsV 
Rs RARs +Rsn) 


(42) 


It can be observed that Eq. (42) is an explicit equation and it has 
only three parameters n, R, and R,,. These parameters can be 
determined by the numerical fitting method. The initial values for 
these three parameters can be obtained as follows: the shunt 
resistance using Eq. (16) and using the linear dependency dV/dI in 
function of (Is-+I—V/Rsn)~'kT/q, the series resistance and the 
ideality factor n of diode can be found by intercept and slope. 
The saturation current and the photogenerated current can be 
calculated using Eqs. (19) and (20). 

Method 11. Garrido-Alzar developed a method to determine 
the solar cell parameters based on the two diode model. The 
method consists of solving the equation system obtained by Eq. (4) 
for four points from the I-V characteristic (Vj, I;), i=1,2,3,4 and 
ng=1, in the zero order approximation for n, and R, [18]. In this 
step, the following parameters can be calculated: Ipn, Ioa, Ior and 
Rsn. To compute n, the following equation is used: 


Voc 
Veg [lor Tpn—loa (ex (Y—1) )—¥e] — In Tor} 


The series resistance is obtained using the following equation: 


Vs +I5Rs V5 +I5Rs 
k= Rata (e aa -1) i (em = 1) me (44) 
sh 


(43) 


The above steps must be repeated until, for example, the value of 
n, is obtained with the desired precision [18]. Here, n,#2 as in 
method 8. 

Method 12. Kiran and Inan developed a method to determine 
the parameters of the solar cell based on the one diode model 
under illumination conditions [19]. They obtained an approxima- 
tion equation from Eq. (3) by eliminating the saturation current 
and considering the shunt resistance as infinite, which is rendered 
as Eq. (45). With such an assumption, it is clear that the method 
provides results which deviate from the expected values when the 
cell is aged. This has a value in the field of study of the PV cell 
aging. 


(45) 


V = Voc—IRs+mV 7 inf ed 


Isc = nVr 


There are only three parameters R, n and Is in Eq. (45). These 
parameters can be computed by fitting the experimental data. The 
photogenerated current can be obtained as lpn ~ Isc. The saturation 
current is given by the following equation: 


= Isc 
exp (4) -1 


Method 13. Sellami and Bouaicha determined the five para- 
meters of the solar cell using the one diode model, the I-V 
characteristic which is measured under illumination conditions 
and for the fitting procedure the genetic algorithm (GA) is used 
[20]. The GA consists of some preliminary decisions: solution 
encoding, evaluation function, initial population generation, selec- 
tion criterion, recombination/reproduction and termination 
criteria. The evaluation function used in GA, derived from Eq. (3) is 


Vi +IRs 
Rsh 


Io (46) 


Vj +IRs 
I(V;, P) =Iph—lo (e ay 1) (47) 
The cost function can be defined in the case of solar cells as the 
sum of the squared difference between the experimental and 
theoretical current values given by Eq. (29). In case of the genetic 
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algorithm, the accuracy criterion being given by the cost function, 
this has to be minimized. The cost function does not involve the 
derivatives and it is determined for each chromosome. The initial 
population of chromosomes IPOP can be calculated as follows [12]: 


IPOP = (hj—I,) random | Nipop, Npar | + lo (48) 


where Npar is the number of parameters, 5 in this case, Nipop 
represents the number of chromosomes, h; and le are the highest 
and the lowest values of the five parameters respectively. The 
selection of the family of chromosomes is made function of the 
cost. The genetic algorithm GA is stopped when the value of the S 
(P) is smaller than the predefined cost minimum. The flow chart 
described by Sellami and Bouaicha is presented in Fig. 3. The same 
method was used by Jervase et al. [21], but they determined the 
parameters of the solar cell for the two diode model. 

Method 14. Ye et al. determined the five or seven parameters of 
the solar cell using the one or two diode model, the I-V char- 
acteristic which is measured under illuminated conditions and the 
particle swarm optimization algorithm (PSO) [22]. The same 
method is described by Qin and Kimball [23]. The PSO algorithm 
allows the determination of the solar cell parameters without 
using the differential equations. This algorithm is an iterative one, 
in which each particle flies in a space with a velocity, which is 
updated as a linear combination by the position, the inertia 
weight, social iteration and self cognition. The velocity and 
position for each particle are updated as follows: 


vi(k+ 1) = w(k)vi(k) +c1r1(k) [p best;(k)—x;(k)] + c2r2(k)[g besti(k)—x;(k)] 
(49) 


xX(kK+1) = vi(k+1)+x;,(k) (50) 


where x; represents the ith particle position (each particle contains 
a vector with the solar cell parameters values), w is the inertia 
weight, v; is the velocity of ith particle, rı and rz are the random 
numbers between 0 and 1, c; and cp are the acceleration constants, 
pbest is the particle’s best known position, gbest represents the 
best known position found by all particles [22,23]. The objective 
function which has to be minimized is 


= piel (51) 


where l, represents the predicted current, Im represents the 
measured current and f; is the fitness function for one sampled 
point. The iteration process is stopped when the fitness function 
reaches minimum conditions. 

Within this concept of solar cell analysis, Sandrolini et al. 
determined the seven parameters of a solar cell using the two 
diode model, the I-V characteristic which is measured under 
illuminated conditions and the PSO algorithm and cluster analysis 
handled together [24]. 

Method 15. Stutenbaeumer and Mesfin used the non-linear 
curve fitting to determine the six parameters of the solar cell using 
the simplified two diode model of a solar cell under dark condi- 
tions [25]: 

V 


I= loa (67-1) Hor (er 1) K (52) 


The simplified two diode model of the solar cell can be obtained 
for small values of R,. In this case, the approximation V>IR, holds. 


Mutate|—» Tetor —>| Stop 
convergence 


Initial Evaluate 
Population cost 


t p J 


Fig. 3. GA flow chart. 
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The I-V characteristic is measured under both illumination and 
dark conditions. In order to reduce the number of fitting para- 
meters the ideality factor of diode ng takes the value 1. The series 
resistance is calculated as follows: 

R; = Vsc—Voc (53) 

Isc 

where Vs. represents the necessary voltage in dark condition to 
obtain the same short circuit current as in illumination conditions. 

The fitting procedure to determine the parameters uses the 
non-linear fitting algorithm developed by Levenberg—Marguardt 
which is implemented in Origin. 

Method 16. Kaminski et al. determined the following para- 
meters of the solar cell Ia, lor, Ng, Nn Rs and Rsn using the two diode 
model in dark conditions [26,27]. The I-V characteristic of the 
solar cell in dark conditions is analyzed on two regions: the first 
represents the higher part where the series resistance and the 
diffusion mechanism are important and the second represents the 
lower part where the shunt resistance and the recombination 
mechanism are important. The first part is used to determine the 
following parameters: loq, Ng and R,. The series resistance and the 
ideality factor of the diode can be determined using the linear 
regression of Eq., (54). 

In(I/h) 


1 
Y= uv et y= l-h `’ x= A (54) 


where (Vj,IJ,) is a point on the first part of the solar cell I-V 
characteristic. If the correlation coefficient is not good, then the 
value of I; must be raised. The saturation current, lọ, can be 
estimated by plotting 1n I in function of V—IR;. The second part 
of the characteristic curve is used to determine the following 
parameters: Ior, Nn, and Rsn. The linear part of the lower part of I-V 
characteristic is fitted by the least square method using Eq. (55). 


I= lor [e —1] (55) 


The shunt resistance is calculated then by Eq. (56). 
-1 
dI lz- I 
Rsn = [O n-Vr a (56) 


Method 17. Arcipiani developed the generalized area method 
[28]. This method allows the determination of the following solar 
cell parameters: the ideality factor of diode, the series and shunt 
resistance. The method consists of calculating the area between 
the OV and OI axes and the I-V characteristic of solar cell, Eq. (57). 
The one diode model was used for the solar cell, using the 
following expression: 


2 
(57) 


where r=R, and g=Rsn. Taking into consideration the following 
approximations gr<1, lo<lpn and lsc%lpn which are generally 
accepted [29], Eq. (57) becomes 


oe {Wei 1 —_IsceVoc—A 
y= Gale (7) vet (Fee (=) Veet se Y= Va 
(58) 


Eq. (58) contains three unknown parameters of the solar cell, r, g 
and n. To determine these parameters Eq. (58) should be rewritten, 
taking I-V at three different illumination levels. Then the resulting 
non-linear system of equations is solved. 

Method 18. The method developed by Jia and Anderson allows 
the determination of the series resistance and the ideality factor of 


I V, 
A= (Iph H Io) (Voc TIsc)t +Isc(1 +81) (vn) +Voeg (Vin—"3*) 
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diode [30]. In this method n is considered variable on the I-V 
characteristic, thus n=1 at Və: and n = 2 at I,,. The one diode 
model is used and the solar cell is illuminated. The series 
resistance can be calculated with Eq. (59) which is obtained using 
the following approximations: lpn ~ Isc, Rsn is infinite and I,->Ip. 


Vin (1/Vr) Usc—Im) {Voce + Vr In [1—(Im/Isc) | }—Im 


R= 
Im (1/Vr) Use—Im){ Voc + Vr ln [1-— (Im/Isc)] }+Im 


(59) 


The value of the ideality factor of the diode at maximum power 
point can be calculated as follows: 


Isc—1 
n= (Vm+InR9/ ( Voc+Vr n ==) (60) 


Method 19. The integral method was developed by Ortiz- 
Conde et al. [31] and Kaminski et al. improved this method [26]. 
The series resistance and the ideality factor of the diode can be 
obtained under dark conditions through the linear regression of 
Eq. (62). This equation is obtained from Eq. (61) through integra- 
tion. The one diode model was used for the solar cell and the shunt 
resistance is considered infinite. 


ik (err —1) (61) 


y=nVr+Rx, X=, y=, | lav (62) 
Vı 

where (Vj,I;) is a point of the solar cell I-V characteristic and I is 

the dark current. The integral can be solved using the means of 

trapezoidal method. 

Method 20. El-Adawil and Al-Nuaim proposed a method to 
determine the series and shut resistances using the one diode 
model [32]. The series resistance can be estimated by using Eq. (3), 
with the approximations shown in Eq. (63) and by choosing two 
points around the knee of I-V characteristic of solar cell (V; Ij), 
i= 1,2, 


IoRene "r >V, Ret <Isc+Io—Hi,i= 1,2, Isc>lo, Isc ™ lpn (63) 


nVr | Isc—l2 Va—V 


Rs = 
$ b-h ls-h bh 
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The shunt resistance can be obtained by the following equa- 
tion: 

Ron + Rs V1—Rs 
RenUse—h)—-hRs nVrh 


(65) 


Thongprona et al. had used the same method before, but for 
both dark and illuminated conditions [33]. 

Method 21. The ideality factor of the diode and the reverse 
current can be determined by using the semi-logarithmic Ise— Voc 
characteristic, as in the equation below. The short circuit current 
and the open circuit voltage must be determined for different 
levels of illumination. In this method, the one diode model is used 
with the following approximations: 


Voc 
Rep oo, Isc © Toh evr >1 (66) 


Under these conditions, Eq. (3) becomes 


1 
In Isc = In Ig tay Vo (67) 


The ideality factor of the diode can be calculated from the slope 
of the semi-logarithmic characteristic and the saturation current is 
calculated using the intercept of the characteristic on Y-axis. 

Priyanka et al. improved this method by taking into considera- 
tion the value of the shunt resistance [34]. Eq. (67) becomes 


Voc 1 
In( Is a) = InI,4 nv; (68) 

Method 22. The method developed by Signal allows the 
determination of the series resistance and the ideality factor of 
diode using the one diode model for one level of illumination [35]. 
The shunt resistance for this method is considered infinite. The 
author developed analytical relations for Pm, Im, Vm and FF in 
function of n and R,. The values for n were chosen 1 or 2. The series 
resistance was considered ideal (R,=0) or real (R, #0). For two 
values of n and Rs, AFF and A(Vm/Voc) were calculated. If n is not 
equal the correct value of n, the value of the series resistance 
obtained from AFF differs from that obtained from A(Vm/Voc). The 
series resistance and the ideality factor of the diode are found by 
varying the value of n until the two values of R, obtained from AFF 
and A(Vm/Voc) are almost equal. 

Method 23. The series and shunt resistance can be determined 
using the one diode model and the I-V characteristic of a solar cell 
in the third and fourth quadrants under illumination conditions. 
This method was developed by Priyanka et al. [34]. The series 
resistance can be calculated using the following equation: 


_iVr In IR- (V-+Vf+I;R) Vy 
l TR I; 


where R=Rs+Rsn, Is, Vp represent I, V in the fourth quadrant and I,, 
V, represent J, V in the third quadrant. The shunt resistance can be 
determined using the Eq. (6) and by making the approximation 
Rs<R, R is determined. The ideality factor of the diode and the 
saturation current must be previously determined using 
method 21. 

The series resistance can be obtained more accurately using 
Eq. (70) by iteration. The first value of the series resistance in 
iterative process is given from Eq. (69). 


Rs 


(69) 


ee a= a 
The values for I, and Ip should be chosen so that 


Method 24. Bowden and Rohatgi developed a method which 
allows the determination of the two parameters of the solar cell, n 
and R, [36]. In this method, two I-V characteristics are used, one 
measured for one sun illumination and the other measured for 
0.1 sun illumination or in the ideal case, at such an illumination 
level that the short circuit current of the solar cell equals the 
difference Isc1ı— 1m1. The series resistance can be calculated using 


_ Voc. = Vint 


Rs (72) 


> Isc, 1=Isc.0,1 
where the pair (Vina, Isc1) is obtained from the I-V characteristic 
measured for one sun illumination and the pair (Voco, Isco1) is 
obtained from the I-V characteristic measured for 0.1 sun illumi- 
nation. The ideality factor of diode can be calculated using 


n= Voc,1 —Voco1 1 
In (Isca /Tsco,1 ) Vr 


Method 25. The series resistance can be determined by the two 
characteristics method which was suggested by Swanson [37], 
described by Wolf and Rauschenbach [38] and reiterated by 
Bashahu and Habyarimana [39]. In this method the I-V character- 
istic of a solar cell is measured at the same temperature for two 


(73) 
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levels of illumination. The series resistance can be calculated using 
the following equation: 


AV _ V2-Vy 


=a e 


(74) 


where the points (Vj, J;), i=1,2 are obtained from the I-V char- 
acteristics, considering the value of AI so that the parallel drawn 
through I, to the Ox axis intersects the characteristic for the 
highest illumination level just a little below the maximum power 
point, see Fig. 4. 

Pysch et al. proposed a modification where at least three I-V 
characteristics were used to determine the series resistance. These 
I-V characteristics can be obtained for three levels of illumination: 
one sun, slightly above and slightly below one sun [40]. 

Method 26. Aratijo and Sanchez developed the area method to 
determine the series resistance of a solar cell [29]. In this method 
the one diode model was used, Eq. (3). Rsp is considered infinite 
and the ideality factor of the diode has the value equal to one. In 
this case, the series resistance can be calculated as follows: 


ToL 6a 


where A represents the area between the I-V characteristic curve 
and the J and V axes. 

Method 27. Cape et al. developed the flash lamp method [41]. 
This method allows the determination of the series resistance 
using Eq. (76). The one diode model was used and Rsn is considered 
infinite. 


R, =R; (Ze) (76) 


where R; represents the load resistance and V; denotes the voltage 
across the load resistance. This method uses two very short pulses 
with duration of 1 ms, to determine V,, and V;. The open circuit 
voltage is determined by using the first light pulse. The voltage 
across the load resistance, which was previously measured with 
precision, is determined using the second light pulse. 

Method 28. The maximum power point method which allows 
the determination of the series resistance was developed by 
Picciano [42,43]. The method uses the one diode model and the 
I-V characteristic is measured for a single illumination level. 
The shut resistance is considered infinite and the ideality factor 
of the diode has a constant value over the entire I-V characteristic. 


I-V Characteristic 


Current [A] 
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Fig. 4. The method with two I-V characteristics for the determination of solar cell 
series resistance. 


The series resistance can be determined as follows: 


Vin 1 
=], eek oo 
where 
B= Um/Use—lm)] + In [Use—lin) /Isc] (78) 


2Vin—Voc 


Method 29. The method developed by Cotfas et al. allows the 
determination of the series resistance using the one diode model 
[44]. This method uses two I-V characteristics - one is measured 
and the other is ideal, Fig. 5. For the ideal case the I-V character- 
istic can be calculated with n=1 or with the real values for n and I, 
obtained with method 21. The series resistance can be calculated 
with Eq. (79). 

R, = AV _ Videa Vm (79) 
Im Im 
where Viaear is obtained on the ideal I-V characteristic for Im. 

Method 30. Agarwal et al. developed a method which allows 
the determination of the series resistance using the one diode 
model and the variation of photogenerated current function of the 
levels of illumination [45]. The shunt resistance for this method is 
considered infinite. The short circuit current can be approximated 
with the photogenerated current for low illumination levels, 
smaller than 750 W/m?. Under low illumination levels, it shows 
linear dependency vs. the illuminations levels. For illumination 
levels higher than 750 W/m? the dependency is under linear. 
The difference between the photogenerated current, by projecting 
the tangent of the Isc at low I, with the Isc curve obtained at normal 
I becomes significant. The series resistance can be found using 
the slope of the semi-logarithmical dependence In(Ipn—Isc) Vs. Isc 
which is linear for illumination levels above 750 W/m: 


Tyn—Isc _ Rslsc 
m( i j-i (80) 


where n and J, can be calculated using method 21 or can be 
approximated with values from ideal cases. 

Method 31. The series resistance of solar cells can be deter- 
mined using the method developed by Polman et al. [46]. This 
method uses the two diode model for the solar cell description. 
The shunt resistance and the ideality factor of diode are taken into 
consideration. Using Eq. (81), the series resistance can be calcu- 
lated as the average between the values obtained when n takes the 
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Fig. 5. The Cotfas method. 
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values 1 and 2. 


-1 
al 1 a 
cad (han) - (zyz) a 


where (dI/oV),, _ v, İS obtained from the slope of the linear fiting 
of a few points around V=V,,. 

Method 32. The method proposed by Cowley and Sze to 
determine the series resistance consists of measuring the semi- 
logarithmical I-V characteristic under dark conditions and using 
the one diode model. The series resistance can be calculated using 
the gap value on the V axis, between the semi logarithmical I-V 
characteristic and the diffusion line [39,47]. 

Method 33. The method allows the determination of the series 
resistance of the solar cells using both measurements under both 
dark and illumination conditions. Aberle et al. [48] proposed the 
determination of the series resistance using Eq. (82): 

Rg = AV_ Vaarkm—Vim (82) 
Im Im 

where Vaarkm represents the voltage measured in dark conditions, 

corresponding to the current Im measured in illumination condi- 

tions. Because the effect of dark series resistance was neglected, 

the error can be larger than 5%. The series resistance is calculated 

with the corrected equation [49] in order to minimize the errors: 


R Vaarkam—Vm—Usc—Im)Rs,dark (83) 
i 
Im 


where Re dark Can be calculated using Eq. (53). 

Method 34. The ideality factor of the diode can be found using 
the analytical method developed by Bayhan and Kavasogl [50]. In 
this case the measurements are performed under dark conditions; 
the authors used the one diode model for the solar cell description 
and the Lambert W-function. The value of n is calculated for two 
regions of I-V characteristic measured under dark conditions. The 
first region is for the voltage interval (0.2, 0.6) and the second 
region is for the voltage interval (0.6, 0.8). 


3. Results and discussion 


This article outlines, discusses and compares a number of 34 
theoretical and experimental approaches in order to determine 
the important dc parameters of solar cells, like Rs, Ryn, n, Ip etc. The 
analysis of the I-V characteristic is an important tool for the 
determination of the above quantities and it can be measured 
under both illumination and darkness conditions. Most methods 
use the I-V characteristic obtained under illumination conditions 
because these are real operating conditions for a solar cell. The 
methods that use the dark measurements are few, for example 
methods 14, 15, 16, 32 and 34. There is, also, a method which uses 
both the dark conditions and the illumination conditions, namely 
method 33 developed by Aberle et al. 

The equation which mathematically describes the behavior of 
solar cells depends on the model chosen for the solar cell. The 
most widely used model is the one diode model because it is the 
simplest and describes sufficiently well the characteristic of most 
solar cells behavior, as it can be seen for all methods except the 
methods 8, 11, 15, 16 and 31. 

The first sixteen methods presented in this paper allow the 
determination of all the parameters: Iph, Io, Rs, Rsn and n. The values 
Of Isc, Voc, Pm and implicitly Vm and Im can be determined from the 
experimental data | 10]. Two or more parameters of a solar cell are 
determined using methods 17-24. The last nine methods allow for 
the determination of one parameter, especially the series resis- 
tance which is the mostly tried parameter for solar cells. 


Most of the methods which allow the determination of all the 
parameters use the least-squares numerical techniques. Eqs. (3) and 
(4) are implicit equations. The non-linear fitting procedure becomes 
complicated. Eq. (3) with the approximation RJ<V becomes an 
explicit equation and the fitting procedure becomes simpler. The 
above approximation holds for a solar cell with small area which 
implies that the current is small or for low level illuminations. Using 
the Lambert W-function, Eq. (3) becomes an explicit equation I= I(V) 
[51] and the fitting procedure simplifies the approach. Moreover, 
using the assumption given by Eq. (15), the explicit equation J=I(V) 
has only three fitting parameters. The extracted process is easier using 
the consecrated software such as Matlab, Mathematica and others. In 
Matlab only a few lines of code are necessary [17]. 

The accuracy of the fitting procedure depends on the following 
factors: the choice of the fitting algorithm, especially on the initial 
values of the fitting parameters and the confidence interval. Devel- 
opment of the generic algorithm and the particle swarm optimization 
may improve the accuracy of the solar cell parameter extraction. A 
high advantage of the method 14 in comparison with the fitting 
methods is that even without very good choice of the values of the 
initial fitting parameters, the final results are very good and the errors 
are very small in comparison to other experimental errors like those 
of methods 25, 26 etc. In method 10 the initial values of the fitting 
parameters are firstly calculated, leading to improvements of the 
results obtained by the fitting process. 

To solve the statistical instability in the fitting algorithm and I- 
V characteristic in concentrator solar cells, Araki and Yamaguchi 
proposed a new mathematical model which describes the solar 
cell behavior: 


VI V-+IRs 
I= lpn I(e 7 1) te (84) 
where e is an error term distributed through a Gaussian pdf 
centered at zero [52]. 

The area method and the generalized area method introduce a 
significant uncertainty brought in through the area limited 
between the two axes and the I-V characteristic. Sometimes, the 
result for Rs comes out negative [39]. 

The methods which use few approximations, such as the 
methods 10, 13 and 14 have to be applied to obtain good values 
for the solar cell parameters. The methods in which the ideality 
factor of diode is given the value equal with 1 and the shunt 
resistance is considered infinite must be used carefully. Method 29 
shows how the value of the series resistance changes when the 
value of n takes the real value or the theoretical value (n=1). 

The solar cell parameters such as: Rs, n, Rsn and I, depend on the 
irradiance levels |17,40,44,53]. Therefore the methods 17, 21 and 
25 which used two or three I-V characteristics obtained for 
different levels of illumination, may be used carefully. The two 
characteristics method 25 can be improved using three or more 
I-V characteristics and R, can be calculated from the inverted slope 
of the linear dependency of the (V;, Ii), i=1,2,3,... 

The series resistance obtained in dark conditions is under- 
estimated. The series resistance is higher due to the larger lateral 
electron flow in the emitter under illuminated conditions, the 
normal operating mode for solar cells. 


4. Conclusions 


In this review paper 34 methods developed over the past 35 
years to determine the main electric parameters of a solar cell 
were critically presented, assessed and discussed. The main para- 
meters in concern are, the reverse diode saturation current, Ip, the 
ideality factor n of diode, the series resistance, R;, the shunt 
resistance, Rsh, and the photocurrent, Ipa. The methods outlined 
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fall into three main groups. The first group includes the methods 
which are based in the theoretical analysis of the I-V curve based 
on techniques which range from simple regression, to non-linear 
regression, stochastic process and fitting based on genetic func- 
tions. This family of approaches may provide the values of all the 
above basic solar cell parameters which play the central role to the 
solar cell performance investigations. The second group includes 
methods which succeed in providing the values of most of the 
important parameters, as the ones above, by dealing with graphi- 
cal analysis of the I-V measured in various illumination conditions, 
either at transient or at steady state. The third group of methods 
uses both sophisticated theoretical analysis of the I-V character- 
istic based on the one or two diode models and graphical analysis 
of the I-V curves. 

These methods are deployed in such a way that they may be 
effectively combined in the most efficient way as to provide the 
values of the above five electric parameters of a solar cell. As result 
of the assessment of the 34 methods outlined in this paper, the 
researcher may choose, according to the experimental set up and 
the tools possessed, some of the inter-related approaches outlined 
above, in order to get a set of values concerning the above solar 
cell parameters useful for further investigation. 
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